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ABSTRACT 



Context. Far-UV photons strongly affect the physical and chemical state of molecular gas in the vicinity of young massive stars. 
Aims. Finding molecular tracers of the presence of FUV radiation fields in the millimeter wavelength domain is desirable, because IR 
diagnostics (PAHs for instance) are not easily accessible towards high extinction line-of-sights. Furthermore, gas phase diagnostics 
provide information on the velocity fields. 

Methods. We have obtained maps of the HCO and H'^CO^ ground state lines towards the Horsehead edge at 5" angular resolution 
with a combination of Plateau de Bure Interferometer (PdBl) and the lRAM-30m telescope observations. These maps have been 
complemented with lRAM-30m observations of several excited transitions at two different positions. 

Results. Bright formyl radical emission delineates the illuminated edge of the nebula, with a faint emission remaining towards the 
shielded molecular core. Viewed from the illuminated star, the HCO emission almost coincides with the PAH and CCH emission. 
HCO reaches a similar abundance than HCO* in the PDR (=: 1-2x10"' with respect to H2). To our knowledge, this is the highest 
HCO abundance ever measured. Pure gas-phase chemistry models fail to reproduce the observed HCO abundance by ~2 orders of 
magnitude, except if reactions of atomic oxygen with carbon radicals abundant in the PDR (i.e., CH2) play a significant role in the 
HCO formation. Alternatively, HCO could be produced in the PDR by non-thermal processes such as photo-processing of ice mantles 
and subsequent photo-desorption of either HCO or H2CO, and further gas phase photodissociation. 

Conclusions. The measured HCO/H'^^CO* abundance ratio is large towards the PDR (-50), and much lower toward the gas shielded 
from FUV radiation (<1). We propose that high HCO abundances (>10"') together with large HCO/H'^CO* abundance ratios (>1) 
are sensitive diagnostics of the presence of active photochemistry induced by FUV radiation. 

Key words. Astrochemistry - ISM clouds - molecules - individual object (Horsehead nebula) - radiative transfer ~ radio lines: ISM 



1. Introduction 

Photodissociation region (PDR) models are used to under- 
, stand the evolution of far-UV (FUV; /iv<13.6eV) illumi- 
' nated matter both in our Galaxy and in external galax- 
ies. Th ese sophis ticated models have been benchmarked re- 
cently ( RoUig et al. 2007) and a r e continuously upgraded (e.g ., 
iGoicoechea & Le Bourlot iGonzalez-Garcfa et al. 20081) . 

Given the large number of physical and chemical processes in- 
cluded in such models, it is necessary to build reference data sets, 
which can be used to test the predictive accuracy of models. Our 
team has contributed to this goal by providing a series of high 
reso lution interferome tric observations of the Horsehead nebula 
(see lPetv et al. 20(j 79, for a summary). Indeed, this source is par- 
ticularly well suited because of its favorabl e orientation and ge - 
ometry, and its moderate distance (~400pc; Habart et al. 20051) . 
We have previously studied the carbon dTevssier et al. 20()4t 
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iPety et al. 2005h and sulfur chemistry (IGoicoechea et al. 20061) . 
and detected the pres ence of a cold de nse core, with active deu- 
terium fractionation dPety et al. 2007ah . 

The formyl radical, HCO, is kno wn to be present i n the inter- 
stellar medium since the late 70's d Snyder et al. 197 6). Snyder 
et al. (1985) give a detailed description of the HCO struc- 
ture and discuss the energy diagram for the lowest energy 
levels. HCO is a bent triatomic asymmetric top with an un- 
paired electron, a-type and fe-type transitions are allowed, with 
a stronger di pole moment (1 .36 Debye) for the a-type tran- 
sitions ( Lands berg et al. 1977h . which are therefore more eas- 
ily detectable. The strongest HCO ground state transitions lie 
at 86.670, 86.708, 86.777 and 86.805 GHz, very close to the 
ground state transition of H'^CO* and to the first rotationally ex- 
cited SiO line. Therefore HCO can be observed simultaneously 
with SiO and H'^CO*. HCO ground state lines have been de- 
tected in the Orion Bar as well as in the den se PDRs NGC 2023, 
NGC 7023 and S 140 dSchilke et al. 20011) . From limited map- 
ping, they have shown that HCO is sharply peaked in the Orion 
Bar PDR, confirming earlier suggestions that HCO is a tracer 
of the cloud illuminated interfaces (de Jong et al. 1980). Garcia- 
Burillo et al. (2002) have mapped HCO' and H'^CO+ in the 
nearby galaxy M82. HCO, CO and the ionized gas present a 
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Table 1. Observation parameters for the maps shown in Fig.[T]and|5](Fig.|5]is available on-line only). The projection center of all 
the maps is azooo = 05''40'"54.27', 52000 = -02°28'00". 



Molecule 


Transition 


Frequency 


Instrument 


Config. 


Beam 


PA 


Vel. Resol. 


Int. Time" 




Noise" 


Obs. date 






GHz 






arcsec 




km s"' 


hours 


K 


K 






1-0 


86.754288 


PdBI 


C&D 


6.76 X 4.65 


13 


0.2 


6.5 


150 


0.10 


2006-2007 


HCO 


lo,i3/2,2-Oo,ol/2,l 


86.670760 


PdBI 


C&D 


6.69 X 4.39 


16 


0.2 


6.5 


150 


0.09 


2006-2007 


CCH 


1,3/2 (2)-0, 1/2(1) 


87.316925 


PdBI 


C&D 


7.24 X 4.99 


54 


0.2 


6.9 


130 


0.07 


2002-2003 



" On-source time computed as if the source were always observed with 6 antennae. * The noise values quoted here are the noises at the mosaic 
phase center (Mosaic noise is inhomogeneous due to primary beam correction; it steeply increases at the mosaic edges). 



Molecule Transition Frequency Instrument # Pix. F^g B^ff Resol. Resol. Int. Time" T,vs Noise Obs. date 

GHz arcsec kms"' hours K mK 



H'^'CO+ 


/= 


1-0 




86.754288 


30m/AB100 


2 


0.95 


0.78 


28.4 


0.2 


2.6/5.0 


133 


69 


2006-2007 


HCO 


lo.l 3/2,2 


-Oo.o 


1/2, 1 


86.670760 


30m/AB100 


2 


0.95 


0.78 


29.9 


0.2 


2.6/5.0 


133 


63 


2006-2007 


HCO 


lo.l 3/2, 1 


-Oo.o 


1/2,0 


86.708360 


30m/AB100 


2 


0.95 


0.78 


29.9 


0.2 


2.6/5.0 


133 


63 


2006-2007 


HCO 


lo.l 1/2, 1 


-Oo.o 


1/2, 1 


86.777460 


30m/AB100 


2 


0.95 


0.78 


29.9 


0.2 


2.6/5.0 


133 


66 


2006-2007 


HCO 


lo.l 1/2,0 


-Oo.o 


1/2, 1 


86.805780 


30m/AB100 


2 


0.95 


0.78 


29.9 


0.2 


2.6/5.0 


133 


66 


2006-2007 



" Two values are given for the integration time: the on-source time and the telescope time. 



Table 2. Observation parameters for the HCO deep integrations shown in Fig.[T] Associated transitions can be found in Table[3] The 
RA and Dec offsets are computed with reference to ffaooo = 05''40"'54.27*, ^2000 = -02°28'00". The positions are also given in the 
coordinate system used to display the maps in Fig.[T]and|5] In this coordinate system, maps are rotated by 14° counter-clockwise 
around the projection center, located at (dx, 5y) - (20", 0"), to bring the illuminated star direction in the horizontal direction and the 
horizontal zero has been set at the PDR edge. 



Position name 


((5RA,(5Dec) 


(Sx,Sy) 




arcsec 


arcsec 


"DCO+ peak" 


(20", 22") 


(44.7", 16.5") 


"HCO peak" 


(-5,0") 


(15.1", 1.2") 



Position 


Frequency 
GHz 


Line area" 
Kkms-' 


Instrument 


F^ff 




Resol. 
arcsec 


Resol. 
km s"' 


Int. Time" 
hours 


K 


Noise 
mK 


Obs. date 


"DCO+ peak" 


86.670760 


0.23 ± 0.009 


30m/B100 


0.95 


0.78 


28.4 


0.27 


0.75/1.5 


134 


11 


2008 




86.708360 


0.12 ±0.009 


30m/B100 


0.95 


0.78 


28.4 


0.27 


0.75/1.5 


134 


11 


2008 


"HCO peak" 


86.670760 


0.52 ± 0.008 


30m/B100 


0.95 


0.78 


28.4 


0.27 


0.75/1.5 


127 


10 


2008 




86.708360 


0.31 ±0.007 


30m/B100 


0.95 


0.78 


28.4 


0.27 


0.75/1.5 


127 


10 


2008 




173.3773770 


0.47 ± 0.023 


30m/C150 


0.93 


0.65 


14.2 


0.067 


2.0/4.0 


667 


66 


2008 




173.4060816 


0.26 ±0.018 


30m/C150 


0.93 


0.65 


14.2 


0.067 


2.0/4.0 


667 


66 


2008 




173.4430648 


0.23 ± 0.020 


30m/C150 


0.93 


0.65 


14.2 


0.067 


2.0/4.0 


667 


66 


2008 




260.0603290 


0.16 ±0.019 


30m/C270 


0.88 


0.46 


9.5 


0.18 


3.0/6.0 


740 


59 


2008 




260.0821920 


0.14 ±0.020 


30m/C270 


0.88 


0.46 


9.5 


0.18 


3.0/6.0 


740 


59 


2008 




260.1335860 


0.12 ±0.017 


30m/C270 


0.88 


0.46 


9.5 


0.18 


3.0/6.0 


740 


59 


2008 




260.1557690 


0.06 ± 0.016 


30m/C270 


0.88 


0.46 


9.5 


0.18 


3.0/6.0 


740 


59 


2008 



" Values obtained from Gaussian fits performed on the spectra using the main beam temperature scale. * Two values are given for the integration 
time: the on-source time and the telescope time. 



nested ring morphology, with the HCO peaks being located fur- 
ther out compared to CO and the ring of H 11 regions. The chem- 
istry of HCO is not well understood. SchiUce et al. (2001) con- 
cluded that it is extremely difficult to understand the observed 
HCO abundance in PDRs with gas phase chemistry alone. As 
a possible way out, they tested the production of HCO by the 
photodissociation of formaldehyde. In this model, H2CO is pro- 
duced in grain mantles, and released by non-thermal photo- 
desorption in the gas phase in the PDR. However, even with 
this favorable hypothesis, the model can not reproduce the 
abundance and spatial distribution of HCO because the photo- 
production is most efficient at an optical depth of a few magni- 
tudes where the photodissociation becomes less effective. 

In this paper, we present maps of the formyl radical ground 
state lines at high angular resolution towards the Horsehead neb- 



ula, and the detection of higher energy level transitions towards 
two particular lines of sights, one in the PDR region and the other 
in the associated dense core. These observations enable us to ac- 
curately study the HCO spatial distribution and abundance. We 
present the observations and data reduction in section |2] while 
the results and HCO abundance are given in section [3] and the 
discussion of HCO chemistry and PDR modeling is given in sec- 
tion |4] 

2. Observations and data reduction 

Tables [T] and |2] summarize the observation parameters for the 
data obtained with the IRAM PdBI and 30m telescopes. The 
HCO ground state lines were observed simultaneously with 
H'^CO^and SiO. Frequency-switched, on-the-fly maps of the 
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H>^CO+ 7=1-0 and HCO ground state lines (see Fig. [S]), ob- 
tained at the IRAM-30m using the A 100 and BlOO 3mm re- 
ceivers (~ 7 mm of water vapor) were used to produce the short- 
spacings needed to complement a 7-field mosaic acquired with 
the 6 PdBI antennae in the CD configuration (baseline lengths 
from 24 to 176 m). The whole PdBI data set will be comprehen- 
sively described in a forthcoming paper studying the fractional 
ionization across the Horsehead edge (Goicoechea et al. 2009, in 
prep.). The CCH data shown in Fig.[T]have been extensively de- 
scribed in Pety et al. 2005. The high resolution HCO loj - Oo o 
data are complemented by observations of the 2o,2 - lo,i and 
3o,3 - 2o,2 multiplets with the IRAM 30m telescope centered on 
the PDR and the dense core. To obtain those deep integration 
spectra, we used the position switching observing mode. The 
on-off cycle duration was 1 minute and the off-position offsets 
were (5RA,6Dec) - (-100", 0"), i-e. the Hn region ionized by 
crOri and free of molecular emission. Position accuracy is esti- 
mated to be about 3" for the 30m data and less than 0.5" for the 
PdBI data. 

The data processing was done with the GILDA^ soft- 
wares (Pety 2005b). The IRAM-30m data were first 
calibrate d to the T*^ scale using the chopper wheel 
method jPenzias & Burrus 1973b . and finally converted to 
main beam temperatures (T^\,) using the forward and main 
beam efficiencies (Feff & Beff) displayed in Table |2] The result- 
ing amplitude accuracy is ~ 10%. Frequency-switched spectra 
were folded using the standard shift-and-add method, after 
baseline subtraction. The resulting spectra were finally gridded 
through convolution by a Gaussian. Position-switched spectra 
were co-added before baseline subtraction. Interferometric 
data and short-spacing data were merged before imaging and 
deconvolution o f the mosaic, u sing standard techniques of 
GILDAS (see e.g. lPetv et al. 2005L for details). 

3. Results and Discussion 

3.1. Spatial distribution 

Fig. [U shows a map of the integrated intensity of the strongest 
HCO line at 86.671 GHz, of the Hi^CO^ 7=1-0 fine and of 
the strongest CCH fine at 87.317 GHz. Fig. |2] displays high 
signal-to-noise ratio spectra of several hyperfine components of 
three HCO rotational transitions towards the HCO and the DCO^ 
emission peaks. 

Most of the formyl radical emission is concentrated in a 
narrow structure, delineating the edge of the Horsehead nebula. 
Low level emission is however detected throughout the nebula, 
including towards the d ense core identifie d by its strong DCO^ 
and H'^CO+ emission dPetv et al. 2007ah . The HCO emission 
is resolved by our PdBI observations. From 2-dimensional 
Gaussian fits of the image, we estimate that the emission width 
is ~13+4" in the plane of the sky. The H'^CO^ emission shows a 
different pattern : most of the signal is associated with the dense 
core behind the photodissociation front, and faint H'^CO^ emis- 
sion detected in the illuminated edge. The CCH emission pattern 
is less extreme than HCO, but shows a similar enhancement in 
the PDR. 

In summary, the morphology of the HCO emission is 
reminiscent of the emission of the PDR tracers, either 
the PAH emission dAbergel et al. 20021) or the small hydro- 
carbons emission, which is strongly enhanced towards the 



' See http://www.iram.fr/IRAMFR/GILDAS for more informa- 
tion about the GILDAS softwares. 



Table 3. Einstein coefficients and upper level energies. 



Molecule 


Transition 


Frequency 








J,F- J', F' 


GHz 


(s ') 


(K) 


HCO 


loi - Ooo 










3/2,2 - 1/2, 1 


86.670760 


4.69 X 10"'' 


4.2 




3/2, 1 - 1/2,0 


86.708360 


4.60 X 10"'' 


4.2 




1/2, 1, 1/2,0 


86.777460 


4.61 X 10"' 


4.2 




1/2,0- 1/2, 1 


86.805780 


4.71 X 10"* 


4.2 




2o2 - loi 










5/2, i - 3/2, 2 


1 1 i.i / /i / lU 


4.51 X lU 


12.5 




5/2,2-3/2, 1 


173.4060816 


4.43 X 10"^ 


12.5 




J/Z, z — i/z, 1 


1 / j.44jUD4o 


i.iy X lU 


iZ.j 




3o3 - 2o2 










7/2,4-5/2,3 


260.0603290 


1.63 X 10"'' 


25.0 




7/2,3-5/2,2 


260.0821920 


1.61 X 10"'' 


25.0 




5/2,3-3/2,2 


260.1335860 


1.45 X 10"'' 


25.0 




5/2,2-3/2, 1 


260.1557690 


1.37 X 10"" 


25.0 


H"CO+ 


7=1-0 


86.754288 


3.2 X 10"' 


4.2 




7=3-2 


260.2553390 


1.3 X 10"3 


25.0 



The line frequenc ies and inte nsities were extracted from the JPL 
jPickett et al. 19 98*) and CDMS jMuller et al. 200lhlMuller et al. 20051) 
molecular spectroscopy data bases for HCO and H'^CO* respectively. 



PDR (tTeyssier et al. 2004t iPetv et al. 20051) . By contrast, the 
HCO emission becomes strikingly faint where the gas is dense 
and shielded from FUV radiation. These reg ions are associate d 
with bright DCO+ and H'^CO^ emission (I Petv et al. 2007al) . 
Our maps therefore confirm that HCO is a PDR species. 

3.2. Column densities and abundances 

3.2.1 . Radiative transfer models of HCO and H^^CO^ 

Einstein coefficients and upper level energies of the studied HCO 
and H'^CO^ lines are given in Table[3] As no collisional cross- 
sections with H2 nor He have been calculated for HCO so far, 
we have computed the HCO column densities assuming a sin- 
gle excitation temperature Tex for all transitions. Nevertheless 
our calculation takes into account thermal, turbulent and opac- 
ity broadening as well as the cosmic m icrowave background 
and line opacity dGoicoechea et al. 2006h . For H^^CO^, detailed 
non-local and non-LTE excitation and radiative transfer calcula- 
tions have been performed using the same approach as in our 
previous PdBI CS and C'^O line analysis (see Appendix in 
Goicoechea et al. 2006). H'^CO^-H2 collisional rate coefficients 
were adapted from those of Flower (1999) for HCO^, and spe- 
cific H'^CO^-electron rates where kindly provided by Faure & 
Tennyson (in prep.). 

3.2.2. Structure of tine PDR in HCO and H^^CO^ 

To get more insight on the spatial variation of the HCO and 
H'^CO^column densities and abundances, we have analyzed a 
cut through the PDR, centered on the "HCO peak" at 6y-Q" (see 
Fig-Ell. The cut clearly shows that HCO is brighter than H'^CO^ 
in the PDR and vice-versa in the dense core. Taking into account 
the different level degeneracies of both transitions (a factor 2.4) 
and the fact that the associated Einstein coefficients Ajj differ by 
a factor ~8 (due to the different permanent dipole moments, see 
TableO, A^(H'^CO+) must be significantly lower than A^(HCO) 
towards the PDR. 
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Fig. 1. High angular resolution maps of the integrated intensity of H'^CO^, HCO, CCH and vibrationally excited H2 emission. 
H'^CO^ and HCO have been observed simultaneously, both with the IRAM-30m and IRAM-PdBI. Maps have been rotated by 
14° counter-clockwise around the projection center, located at (6x, 6y) - (20", 0"), to bring the illuminated star direction in the 
horizontal direction and the horizontal zero has been set at the PDR edge. The emission of all lines is integrated between 10.1 and 
11.1 km s Displayed integrated intensities are expressed in the main beam temperature scale. Contour levels are displayed on the 
grey scale lookup tables. The red vertical line shows the PDR edge and the green crosses show the positions (DC O^ and HCO peaks ) 
where deep integrations have been performed at IRAM-30m (see Fig.|2]i. The H2 map is taken from (author?) dHabart et al. 2005h . 



We modeled the PDR as an edge-on cloud inclined by ~5° 
relative to the line-of-sight. We have chosen a cloud depth of 
~0.1pc, which implies an extinction of Ai/^20 mag for the 
considered densities towards the "HCO peak". These parame- 
ters are the best geometrical d escription of the Horsehead PDR- 
edge (e.g.. iHabart et al. 2005h and also reproduce the observed 
1.2 mm continuum emission intensity. The details of this mod- 
eling will be presented in Goicoechea et al. (2009). In the fol- 
lowing, we describe in details the determination of the column 
densities and abundances for two particular positions, namely 
the "HCO peak" and the "DCO^ peak" (offsets relative to the 
map center can be found in Table |2]). 

3.2.3. HCO column densities 

We used the three detected rotational transitions of HCO (each 
with several hyperfine components, see Fig. |2]l to estimate the 
HCO column densities in the direction of the "HCO" peak. We 
have taken into account the varying beam dilution factors of the 
HCO emission at the "HCO peak" by modeling the HCO emis- 
sion as a Gaussian filament of ~12" width in the 6x direction, 
and infinite in the by direction. The filling factors at 260, 173 
and 87 GHz are thus ~0.8, 0.6 and 0.4, respectively. 

A satisfactory fit of the IRAM-30m data towards the "HCO 
peak" is obtained for Tex -5 K and a turbulent velocity dis- 
persion of cr=0.225kms-' (FWHM= 2.355 x cr). Line profiles 
are reproduced for A^(HCO) - 3.2 x 10'^ cm"^ (see red solid 
curves in Fig. |2|. The most intense HCO lines at 86.67 and 
173.38 GHz become marginally optically thick at this column 
density (r > 1). Therefore, opacity corrections need to be taken 
into account. We checked that the low value of Tex (subthermal 
excitation as -60 K) is consistent with detailed excitation cal- 
culations carried for H'-'CO^ in the PDR which are described 
below. 

Because the HCO signals are weaker towards the 
"DCO^peak", we only detected 2 hyperfine components of the 
loi - Ooo transition. Assuming extended emission and the same 
excitation temperature as for the "HCO peak", 5 K, we fit the ob- 
served lines with a column density of 4.6 x 10'^ cm"^ (blue solid 
lines in Fig. |2]i. Both HCO lines are optically thin at this posi- 
tion. This simple analysis shows that the HCO column density 



Table 4. Inferred column densities and abundances with respect 
to molecular hydrogen, e.g. - N(X)/N(H2)- 



Molecule 


Method 


HCO peak 


DCO+ peak 


A'(H2)[cm-^] 


1.2 mm cont. 


1.9 X 10"" 


2.9 X 10"" 


A'(HCO)[cm--] 


T-ex = 5 K 


3.2 X 10' ' 


4.6 X 10'" 


A'(H'3C0+)[cm--] 


Full excitation 


5.8 X 10" 


5.0 X 10'" * 


A'(HC0+)[cm-2] 


'2C/'3C=60 


3.5 X 10'3 


3.0 X 10'" 


A^(HCO) 




1.7 X 10-'' 


1.6 X 10-'" T 


;r(H"co+) 




3.1 X 10-" 


1.7 X 10-'" 


;r(HCO+) 




1.8 X 10"' 


1.0 X lO-** 



• Pety et al. (2007a) ' 1.7 x 10-' if HCO arises only from the cloud 
surface (A^ - 3). 



is ~ 7 times larger at the "HCO peak" in the PDR, than towards 
the dense cold core. 

3.2.4. H'^CO+ column densities 

Both the H'^CO+ 7=3-2 and 1-0 line profiles at the "HCO peak" 
are fitted with n(H2) ^5x10"* cm-^, Tk ^60K and e /H ^5x10-^ 
(as predicted by the PDR models below). The required column 
density is A^(H'^CO^)=5.8xI0'' cm"^. For those conditions, the 
excitation temperature, Tgx, of the 7=3-2 transition varies from 
^4 to 6 K, which supports the single-Tex models of HCO. Both 
H'^CO^ lines are optically thin towards the "HCO peak". 

The H'^CO^ line emission towards the "DCO^ peak" has 
been studied by Pety et al. (2007a). Both H'^CO^ Unes are mod- 
erately optically thick towards the core, and the H'^CO^ column 
density is A^(H'''CO^)- 5.0x 10'~ cm-^, which represents an en- 
hancement of nearly one order of magnitude relative to the PDR. 
According to our 1 .2 mm continuum map, the extinction towards 
the core is Ay >30 mag compared to 20 mag in the PDR. The 
H'^CO^column density enhancement therefore corresponds to a 
true abundance enhancement. 

3.2.5. Comparison of HCO and H'^GO+ abundances 

Table|4]summarizes the inferred HCO and H'^CO^ column den- 
sities and abundances towards the 2 selected positions : the 
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5 



HCO 



HCO 



HCO 



0.8 : 

PDR 
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Fig. 2. IRAM-30m observations (histograms) of several HCO hyperfine components of t he lni-Onn, 2n?-ln i and 3o3-2o2 rotational 
transitions towards the PDR ("HCO peak") and towards the dense core ("DCO+ peak") dPetv et al. 2007ah . Solid lines are sinsle- 
Tex radiative transfer models of the FDR-filament (red curves) and line-of-sight cloud surface (blue curves). A sketch of the HCO 
rotational energy levels is also shown (right comer). 



"HCO peak" in the PDR and the "DCO+ peak" in the FUV- 
shielded core. Both species exhibit strong variations of their 
column densities and abundances relative to H2 between the 
PDR and the shielded region. In the PDR, we found that both 
the HCO abundance relative to H2 Ct-(HCO)-l-2 xlO '*) and 
the HCO/H'''CO+ column density ratio (^50) are high. These 
figures are higher than all previously published measurements 
(at lower angular resolution). Besides, the formyl radical and 
HCO^ reach similar abundances in the PDR. 

The situation is reversed towards the "DCO^ peak", i.e. the 
observed HCO/H'^'CO^ column density ratio is lower (wl) than 
towards the "HCO peak" . Nevertheless, while the bulk of the 
observed H'^CO^ emission arises from cold and shielded gas, 
the origin of HCO emission is less clear HCO could either (/) 
coexist with H'^CO^ or (//) arise predominantly from the line- 
of-sight cloud surface. In the former case, our observations show 
that the HCO abundance drops by one order of magnitude be- 
tween the PDR and the dense core environment. However, it is 
possible that the abundance variation is even more pronounced, 
if the detected HCO emission arises from the line of sight cloud 
surface. We have estimated the depth of the cloud layer, assum- 
ing that HCO keeps the "PDR abundance" in this foreground 
layer : a cloud surface layer of Av ^3 (illuminated by the mean 
FUV radiation field around the region) also reproduces the ob- 
served HCO lines towards the cold dense core (blue solid lines 
in Fig.©. 

In this case, both the HCO abundance and the HCO/H'^CO+ 
abundance ratio in the dense core itself will be even lower than 
listed in table H] We have tried to discriminate both scenar- 
ios by comparing the HCO loi - Ooo (7=3/2-1/2, F=2-l) and 
H'^CO^ 7=1-0 line profiles towards this position. Both lines 
have been observed simultaneously with the IRAM-30m tele- 
scope. Because of their very similar frequencies (~86.7 GHz), 
the beam profile and angular resolution is eff'ectively the same. 




1" 12 ^LSH [kms ] 



Fig. 3. Observations along a horizontal cut through "the HCO 
peak" (histograms). The H'^CO+ 7=1-0 and HCO loi - Ooo 
lines were mapped with the PdBI at an angular resolution of 
6.8", whereas the H'^CO"" 7=3-2 fine was mapped with HERA- 
30m (and smoothed to a spatial resolution of 13.5"). Radiative 
transfer models of an edge-on cloud with a line of sight ex- 
tinction of Av-2Q, inclined 5° relative to the line of sight for 
HCO (red curve), and H'^CO^ (blue curves) are shown. The 
single-Tex HCO model assumes a 12" width filament with a 
column density of 3.2 x lO'-' cm"^, while A^(HCO) is 4.6 x 
10'^ cm"^ behind the filament. The H'^CO^ model assumes 
a constant density of n(H2)=5xl0^cm"^ with 7'k=60K and 
A^(Hi3CO+)=5.8 X 10" cm-2 for ^x<35"; and rk=10K and 
A^(Hi3CO+)=7.6xlO" cm-^foT6x>35". Modeled line profiles 
have been convolved with an appropriate Gaussian beam cor- 
responding to each PdBI synthesized beam or 30m main beam 
resolution. 



In this situation, any difference in the measured linewidths re- 
flects real differences in the gas kinematics and turbulence of 
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Fig. 4. Photochemical models of a unidimensional PDR. Upper panels show the density gradient («// = n(H) + 2n(H2) in cm"^) used 
in the calculation. Middle panels show the predicted HCO and H'^CO^ abundances (relative to «//)■ The H'-^CO^ abundance inferred 
from observations in the cold core ("the DCO^ peak", see the offsets in Table |2]i is shown with blue lines. The HCO abundance 
inferred from observations in the PDR ("the HCO peak", see the offsets in Table|2|i is shown with red lines. Lower panels show the 
HCO/H'^CO^ abundance ratio predicted by the models whereas the HCO/H'^'CO^ column density ratio inferred from observations 
is shown as blue arrows and red lines (for the cold core and PDR respectively). Every panel compares two different models: Left- 
side models show a standard chemistry (dashed curves) versus the same network upgraded with the addition of the H2CO + photon 
— > HCO + H photodissociation (solid curves). Riglit-side models show the previous upgraded standard model (solid curves) versus 
a chemistry that adds the O + CH2 — > HCO + H reaction with a rate of 5.01x10 " cm^ s ' (dotted curves). The inclusion of the O 
+ CH2 reaction has almost no effect on H'^CO^ for the physical conditions prevailing in the Horsehead, but triggers an increases of 
the HCO abundance in the PDR by two orders of magnitude. 



the regions where the line profiles are formed. Gaussian fits of 
the HCO and H'^CO+ lines towards "the DCO+ peak" provides 
line widths of Av(HCO) = 0.8 1+0.06 km s ' and Av(Hi3CO+) 
= 0.60+0.01 kms-i. Therefore, even if the H'^CO^ 7=1-0 line 
are slightly broadened by opacity and do not represent the true 
line of sight velocity dispersion, HCO lines are broader at the 
3cr level of confidence. This remarkable difference supports sce- 
nario (ii) where the H'^CO^ line emission towards the "the 
DCO^ peak" arises from the quiescent, cold and dense core, 
whereas HCO, in the same line of sight, arises predominantly 
from the warmer and more turbulent outer cloud layers. We note 
that the presence of a foreground layer of more diffuse mate- 
rial (Ay ~2 mag) was already introduced by Goicoechea et al. 
(2006), to fit the CS J -2-1 scattered line emission. The anal- 
ysis of CO 7=4-3 and CI Pq maps led also (author?) 

dPhilipp et al. 2006 ) to propose the presence of a diffuse enve- 
lope, with Ay ~2 mag, and which contributes to roughly the 
about half the mass of the dense filament traced by C'^O and 
the dust continuum emission. The hypothesis of a surface layer 
of HCO is therefore consistent with previous modeling of molec- 
ular emission of the horsehead. 

We conclude 1) that HCO and HCO^ have similar abun- 
dances in the PDR, and 2) that the HCO abundance drops by at 
least one order of magnitude between the dense and warm PDR 
region and the cold and shielded DCO^ core. 



4. HCO chemistry 

4.1. Gas-phase formation: PDR models 

In order to understand the HCO and H'^'CO^ abundances and 
HCO/H'^CO^ column density ratio infeiTed from observations, 
we have modeled the steady state gas phase chemistry in the 
Horsehead edge. The density distribution in the PDR is well 
represented by a density gradient nniSx) oc where 6x is the 
distance from the edge towards the cloud interior and nn = 
n(H) + 2n(H2) (see the top panels of Fig. |4]i. The density 
reaches a cons tant nn value of 2x10^ cm~^ in an equi valent 
length of -10" (iHabart et al. 2005HGoicoechea et al. 2006 ). The 
cloud edge is illuminated by a FUV field 60 times the mean 
interstellar radiation field (Go = 60 in Draine units). We used 
the Meudon PDR cod^E a phot ochemica l model of a unidi- 
mensional PDR (seelLe Bourlot et al. IQQBtFLe Petit et al. 20061; 
iGoicoechea & Le Bourlot 20 07. for a detailed description). Our 
standard chemical network is based on a modified version of 
Ohio State University (os u) gas-phase network, updated for 
photochemical studies (see ^Goicoechea et al. 2006|). I t also in- 
cludes '^C fractionation reactions (Graedel et aL and spe- 
cific computation of the '^CO photodissociation rate as a func- 
tion of depth. The ionization rate due to cosmic rays in the 



- Publicly available at http://aristote.obspm.fr/MIS/ 
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models is ^=5x10 '^s Following our previous works, we 
chose the following elemental gas phase abundances: He/H=0. 1, 
O/H=3xl0-^ C/H=1.4xl0-^ N/H=8xlO-^ S/H=3.5xlO-<', 
^3c/H=2.3xlO-^ Si/H= 1.7x10-** and Fe/H^LOxlO"''. 

In Fig.m we investigate the main gas-phase formation routes 
for HCO in a series of models "testing" different pathways lead- 
ing to the formation of HCO. HCO and H'^CO+ predictions 
are shown in Figure |4] (middle panels). As a first result, note 
that in all models the HCO abundance peaks near the cloud 
surface at Ay ^1.5 (5x^14") where the ionization fraction is 
high (e"/H~5xlO"^). Due to the low abundance of metals in the 
model (as represented by the low abundance of Fe), the ioniza- 
tion fraction in the shielded regions is low (e /H <10 **), and 
therefore the H'-'CO^ predictions matches the observed values 
(Goicoechea et al. 2009). Besides, a low metalicity reduces the 
efficiency of charge exchange reactions of HCO^ with metals, 
e.g.. 



Fe + HCO"^ ^ HCO + Fe+ 



which are the main gas-phase formation route of HCO in the 
FUV-shielded gas in our models. Hence, the HCO abundance re- 
mains low inside the core. Nevertheless, despite that such mod- 
els do reproduce the observed HCO distribution, which clearly 
peaks at the PDR position, the predicted absolute HCO abun- 
dances can vary by orders of magnitude depending of the domi- 
nant formation route. 

In our standard model (left-side models : dashed curves), 
the formation of HCO in the PDR is dominated by the disso- 
ciative recombination of H2CO^, while its destruction is domi- 
nated by photodissociation. Even if the predicted HCO/H'^CO^ 
abundance ratio satisfactorily reproduces the value inferred from 
observations, the predicted HCO abundance peak is ~3 orders 
of magnitude lower than observed. In order to increase the gas- 
phase formation of the HCO in the PDR we have added a new 
channel in the photodissociation of formaldehyde, the produc- 
tion HCO, in addition to the normal channel producing CO : 



H2CO + photon HCO + H 



(2) 



This channel is generally not included in standard chemical net- 
works but very likely exists (Troe2007: Yin et al. 2007). We 
included this process with an unattenuated photodissociation 
rate of /cj,i5(H2CO)=10 s ' and a depth dependence given by 
exp(-1.74 Ay). This is the same rate as the one given by van 
Dishoeck (1988) for the photodissociation of H2CO producing 
CO, which is explicitly calculated for the Draine (1978) radia- 
tion field. Model results are shown in Fig. |4] (left-side models: 
solid curves). The inclusion of Reaction |2l which becomes the 
dominant HCO formation route, increases the HCO abundance 
in the PDR by a factor ~5. But the HCO production rate is still 
too low to reproduce the abundance determined from observa- 
tions. 

Another plausible possibility to increase the HCO abundance 
in the PDR by pure gas-phase processes is to include additional 
reactions of atomic oxygen with carbon radicals that reach high 
abundances only in the PDR. Among the investigated reactions, 
the most critical one. 



O -I- CH2 ^ HCO + H 



(3) 



is known to proceed with a relatively fast rate at high tem- 
peratures (5.01x10"" cm^s-i at =1200- 1800 K; Tsuboi & 
Hashimoto 1981). This is the rate recommended by N IST 
dMallard et al. 1994 and UMIST2006 dWoodall et al. 200% and 
that we adopt for our lower temperature domain (~10-200K). 



Model predictions are shown in Fig.|4](right-side models: dotted 
curves). Whereas the predicted HCO abundance in the shielded 
gas remains almost the same, the HCO abundance is dramati- 
cally increased in the PDR (by a factor of ~125) and the O + 
CH2 reaction becomes the HCO dominant production reaction. 
Therefore, such a pure gas-phase model adding reactions |2] and 
[3] not only reproduces the H'^CO^ abundance in shielded core, 
but also reproduces the observed HCO absolute abundances in 
the PDR. In this picture, the enhanced HCO abundance that we 
observe in the Horsehead PDR edge would be fully determined 
by the gas-phase chemical path : 



ch: 



CH+ 



^ CH2 ^ HCO 



(4) 



The validity of the rate of Reaction[3]used in our PDR model 
remains, of course, to be confirmed theoretically or experimen- 
tally at the typical ISM temperatures (10 to 200 K). 



(1) 4.2. Other routes for HCO formation: Grain photodesorption 



If Reaction [3] is not included in the chemical network, the pre- 
dicted HCO abundance is ~2 orders of magnitude below the ob- 
served value towards the PDR. As a consequence, the presence 
of HCO in the gas-phase should be linked to grain mantles for- 
mation routes, and subsequent desorption processes (not taken 
into account in our modeling). In particular, SchiUce et al. (2001) 
proposed that HCO could result from H2CO photodissociation, 
if large quantities of formaldehyde are formed on grain man- 
tles and then released in the gas phase. Even with this assump- 
tion, their model could not reproduce the observed HCO abun- 
dance in highly illuminated PDRs such as the Orion Bar The 
weaker FUV-radiation field in the Horsehead, but large density, 
prevent dust grains to acquire high temperatures over large spa- 
tial scales. In fact, both gas and grains cool down below ~30 K 
in ^10"-20" (or Ay ^1-2) as the FUV-radiation field is attenu- 
ated. Therefore, thermal desorption of dust ice-mantles (presum- 
ably formed before cr-Orionis ignited and started to illuminate 
the nebula) should play a negligible role. Hence a non-thermal 
desorption mechanism should be considered to produce the high 
abundance of HCO observed in the gas phase. This mechanism 
could either produce HCO directly or a precursor molecule such 
as formaldehyde. 

Since high HCO abundances are only observed in the 
PDR, FUV induced ice-mantle photo -desorption (with rates 
that roughly scales with the FU V-radiation field strength) 



seems the best candidate (e.g., IWillacv & Williams 1993 



iBergin et al. 19951) . Laboratory experiments have shown that 
HCO radicals are produced in irradiated, methanol con- 



tain ing, ice mantles (Bernstein et al. 19951 : iMoore et al. 20011: 
Bennett & Kaiser 2007h . The formyl radical could be formed 
through the hydrogenation of CO in the solid phase. It is an 
important intermediate radical in the synthesis of more com- 
plex o rganic molecules such as methyl formate or glycolalde- 
hyde feennett & Kaiser 2007h . However, the efficiency of the 
production of radicals in FUV uradiated ices remains uncer- 
tain, and very likely depends on the ice-mantle composition. 
The formation of species like formaldehyde and methanol in 
CO-ice exposed to H-atom bombardment has be en reported 
by different groups jHiraoka et al. 19941: IWatanabe et al. 20021 : 
iLinnartz et al. 20071) . further confirming the importance of HCO 
as an intermediate product in the synthesis of organic molecules 
in ices. Indeed, hydrogenation reactions of CO-ice, which 
form HCO, H2CO, CH 3 O and CH3OH in gr ain mantles (e.g., 
iTielens & Whittet 19971 ICharnlev et al. 19971) . are one impor- 
tant path which warrants further studies. 
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To compare with our observations, we further need to un- 
derstand how the radicals are released in the gas phase, either 
directly during the photo-processing, or following FUV induced 
photo-desorption. Recent laboratory measurements have started 
to shed light on the efficiency of photo-desorption, which de- 
pends on the ice composition and molecule to be desorbed. 
For species such as CO, the rate of photo-desorbed molecules 
per FUV photon is much larger than previously thought (e.g., 
ibberg et al. 20071) . Similar experiments are required to constrain 
the formation rate of the various species that can form in inter- 
stellar ices and to determine their photo-desorption rates. 

We can use the measured gas phase abundance of HCO to 
constrain the efficiency of photo-desorption. We assume that the 
PDR is at steady state, and that the main HCO formation mecha- 
nism is non thermal photo-desorption from grain mantles (with a 
Fhco rate), while the main destruction mechanism is gas-phase 
photodissociation (with a Dhco rate), therefore : 



Dhco = Go ^diss(HCO);r(HCO) n(H2) [cm-^ s"' ] 



(5) 



Fhco = Go ^pd(HCO);ir(HCOice) "^"/^'''^ n(H2)[cm-3 s"'] (6) 

n(H2) 

where xiUCO) is the gas phase abundance of HCO relative to 
H2,;if(HCOice) is the solid phase abundance relative to water ice, 
and n(H20ice)/n(H2) is the fraction of water in the solid phase 
relative to the total gas density. /(•diss(HCO) and /Cpd(HCO) are the 
HCO photodissociation and photo-desorption rates respectively. 
By equaling the formation and destruction rates, we get : 



^pd(HCO) = ^diss(HCO) 



;r(HCO) «(H2) 



;kr(HCOiee) «(H20ice) 



or 



KpaiHCO) 



10 



,2 ^diss(HCO) ;r(HCO)/10-'^ IQ-UUz) 

10-9 ;r(HCOice)/10-2 n(H20ice) 



(7) 



(8) 



where we have used typical figures for the HCO abundance in 
the gas phase (-lO"'^, see above) and solid phase (-10"^ see e.g. 
Bennet & Kaiser 2007) and for the amount of oxygen present as 
water ice in grain mantles. 

Assuming standard ISM grains with a radius of 0. 1 fim the 
required photodesorption efficiency (or yield) Fprf(HCO): 



^pd(HCO) 



Go exp(—2Av) na^ 



[molecules photon 



(9) 



(see e.g., Id'Hendecourtet al. 1985l:lBerginetal. 1995h converts 
to Yp^qAQQi) ^ \^-^ molecules per photon (for the FUV radi- 
ation field in the Horsehead and Ay ^1.5, where HCO peaks). 
Therefore, the production of HCO in the gas phase from photo- 
desorption of formyl radicals could be a valid alternative to gas 
phase production, if the photo-desorption efficiency is high and 
HCO abundant in the ice mantles. This mechanism also requires 
further laboratory and theoretical studies. 

Because the formyl radical is closely related to formalde- 
hyde and methanol and the three species are likely to coexist in 
the ice mantles, a combined analysis of the H2CO, CH3OH and 
HCO line emissions towards the Horsehead nebula (PDR and 
cores) is needed to provide more information on the relative effi- 
ciencies of gas-phase and solid-phase routes in the formation of 
complex organic molecules in environments dominated by FUV- 
radiation. This will be the subject of a future paper. 



5. Summary and conclusions 

We have presented interferometric and single-dish data showing 
the spatial distribution of the formyl radical rotational lines in 
the Horsehead PDR and associated dense core. The HCO emis- 
sion delineates the illuminated edge of the nebula and coincides 
with the PAH and hydrocarbon emission. HCO and HCO^ reach 
similar abundances (^1-2x10 9) in these PDR regions where 
the chemistry is dominated by the presence of FUV photons. 
For the physical conditions prevailing in the Horsehead edge, 
pure gas-phase chemistry is able to reproduce the observed HCO 
abundances (high in the PDR, low in the shielded core) if the 
O -I- CH2 — > HCO + H reaction is included in the models. This 
reaction connects the high abundance of HCO, through its for- 
mation from carbon radicals, with the availability of C^ in the 
PDR. 

The different linewidths of HCO and H'^CO+ in the line of 
sight towards the "DCO^ peak" suggest that the H'^CO^ line 
emission arises from the quiescent, cold and dense gas com- 
pletely shielded from the FUV radiation, whereas HCO predom- 
inantly arises from the outer surface of the cloud (its illuminated 
skin). As a result we propose the HCO/H'^CO^ abundance ra- 
tio, and the HCO abundance itself (if >10 '"), as sensitive di- 
agnostics of the presence of FUV radiation fields. In particular, 
regions where the HCO/H'^CO^ abundance ratio (or intensity 
ratio if lines are optically thin) is greater than ^1 should reflect 
ongoing FUV-photochemistry. 

Given the rich HCO spectrum and the possibility to map its 
bright millimeter line emission with interferometers, we propose 
HCO-H2 as a very interesting molecular system for calculations 
of the ah initio inelastic collisional rates. 
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Fig. 5. Medium angular resolution maps of the integrated intensity of the 4 hyperfine components of the fundamental transition of 
HCO. These lines have been observed simultaneously at IRAM-30m. Maps have been rotated by 14° counter-clockwise around the 
projection center, located at (6x, 6y) - (20", 0"), to bring the illuminated star direction in the horizontal direction and the horizontal 
zero has been set at the PDR edge. The emission of all lines is integrated between 9.6 and 11.4 kms '. Displayed integrated 
intensities are expressed in the main beam temperature scale. Contour levels are displayed on the grey scale lookup tables. The red 
vertical line shows the PDR edge and the green crosses show the positions (DCO^ and HCO peaks) where deep integrations have 
been performed at IRAM-30m (see Fig.|2]i. 



